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ABSTRACT 

Collisional ring galaxies are the outcome of nearly axisymmetric high-speed encounters 
between a disc and an intruder galaxy. We investigate the properties of collisional 
ring galaxies as a function of the impact parameter, the initial relative velocity and 
the inclination angle. We employ new adaptive mesh refinement simulations to trace 
the evolution with time of both stars and gas, taking into account star formation and 
supernova feedback. Axisymmetric encounters produce circular primary rings followed 
by smaller secondary rings, while off-centre interactions produce asymmetric rings with 
displaced nuclei. We propose an analytical treatment of the disc warping induced by 
an inclination angle greater then zero. The star formation history of our models is 
mainly influenced by the impact parameter: axisymmetric collisions induce impulsive 
short-lived starburst episodes, whereas off-centre encounters produce long-lived star 
formation. We compute synthetic colour maps of our models and we find that rings 
have aB-F colour typically ~ 0.2 mag bluer than the inner and outer disc, in 
agreement with observations. 

Key words: galaxies: interactions - galaxies: peculiar - galaxies: star formation - 
methods: iV-body simulations. 



1 INTRODUCTION 

Collisional ring galaxies (CRGs) consist of 
or more sharply defined star-forming rings, 
erally (but not alway s) wit h a nucleus 
l|Madore. Nelson fe PetrTliol 120091 ). Theoretical 



els, both analytica l (e.g . Struck-Marcelll 



gen- 
inside 
mod- 



Applcton & Struck-Marcel 



(e.g. iLvnds fc Toomrel 11976: iThevs fc Spiegell 



Struck-Marcell fc Lotan 199C ; Gerber fc Lambl 



1996) 



1990 



1994 



and numerical 



Appleton fc Struck-Marcelll " 19871: Appleton fc James! 



Hernquist fc Weill 1 19931: IStruck-Marcell fc Higdonl Il995 
Gerber. Lamb fc Balsaral Il994l. Il99rj; iHorellou fc Combe* 



1977 



1990 



bed 
>clli 



20011 : ID'Onghia, Mapelli fc Moord [20081: iGhosh fc Mapellil 
20081: IMapelli et al.1 l2008al . l2008bl : iMapelli fc Maverl 120121 : 
Smith et al.1 120121 ). suggest that CRGs are the result of 



nearly axisymmetric high-speed encounters between a disk 
galaxy (the 'target' or 'primary' galaxy) and an intruder 
galaxy. After the interaction, radially expanding density 
waves form because of the crowding of star orbits in the 
target disk. This causes gas compression and the triggering 
of a starbur s t episode along the e x panding ring (e.g . 



ex pa 

iHigdonl 1 19951 : iHigdon fc Wallinl 1 19971 : iMavva et alTbOOSl : 
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Bizyacv, Mois eev fc Vorobvov 2007; R appaport et al. 2010; 



Fogartv et al.ll201ll ) 



CRGs are unique laboratories for the study of galaxy 
collisions because of their simple interaction geometry. 
For example, the kinematics of the ring provides in- 
formation about the dyna mi cs of the interactio n (e .g. 



Fosburv fc Hawardenl 19771 : 
Charmandaris fc Appleton! 



.Few. Arp fc Madorel 
Il994l : iBizvaev et ail 



1982; 



2007; 



Fogartv et al.l 201ll ). This information is used to constrain 



Af-body/hydrodynamical models reproducing the formation 



1 19931 : 

and 



iHernauist & Weill 19931: 


Struck-Marcell & Higdon 


Horellou & Combesl 12001 


: Mapelli et al. 2008a) 


Arp 147 (e.g. iGerber. Lamb & Balsaral (1992: 



IMapelli fc Maverll201S 

Numerical simulations are helpful tools to explore the 
parameter space of the collision and to disenta ngle the effect 
of eac h parameter on the evolution of CRGs. IGerber et al.l 
( 1996) study the properties of CRGs formed in axisymmet- 
ric encounters with different intruder-to-target mass ratios. 
They show that, as the intruder mass increases, the ring 
becomes sharper, exp ands more rapidly and th e secondary 
ring is less developed. IMapelli fc Maverl (|2012l ) explore the 
connection between gas fraction and star formation (SF) in 
simulations of empty CRGs (i.e. CRGs without a central 
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nucleus, such as Arp 147), and find that the gas fraction 
strong ly influences t he pe ak of SF rate (SFR) after the star- 
burst. ISmith et al.l (|2012l ) perform a large parameter study 
with the aim of reproducing the d ynamical properties of the 
Auriga's Wheel |Conn et al.ll20lil) . but they do not take into 
account the effects of the interaction parameters on the SF. 

In this paper, we investigate how the geometry of the 
impact affects the morphology and SF history of CRGs. We 
employ a suite of adaptive mesh refinement (AMR) numer- 
ical simulations to study systematically the properties of 
CRGs as a function of the interaction parameters that de- 
fine the geometry of the encounter. This is the first time 
that the AMR technique is applied to a wide grid of simu- 
lations studying the formation of CRGs. Our choice is par- 
ticularly important, as the AMR technique ensures a better 
treatment of hydrodynamical instabiliti es and shocks with 
respect to particle-based schemes (e.g. lAgertz et al.l 120071 ; 
iPricd 12008). The paper is organized as follows: in Section 
[5] we describe the main characteristics of our suite of sim- 
ulations and briefly discuss the parameters of interest. In 
Section[3]we present the results, and focus on the morphol- 
ogy and SF history of our models. We also extract synthetic 
photometric information about the colour of the ring and we 
compare it with observations. In SectionUwe summarize our 
main findings. 



2 NUMERICAL SIMULATIONS 
2.1 Numerical code 

Our simulations of galaxy interactio ns were perform ed with 
the public AMR code RAMSES (|Tevssierj 120021 '). RAM- 
SES solves the fluid equations on an adaptive Cartesian 
grid by means of a second-order Godunov method cou- 
pled with a Riemann solver. This method ensures a better 
treatment of hydrodynamical instabiliti es and shocks with 
respect to particle-based schemes (e.g. lAgertz et al.l 120071 ; 
|Pricdl2008l ). Gravitational interactions are calculated using 
Particle Mesh (PM) techniques for iV-body collisionless com- 
ponents (stars and dark matter particles), solving the Pois- 
so n equation on the AMR grid with the multigrid strategy 
of iGuillet fe Tevssierl (|201lh . 

All the simulations presented in this paper adopt a cu- 
bic box of 960 kpc per edge, with a static coarse grid com- 
posed of 64 3 cells. We used 7 additional levels of refinement 
to reach the maximum resolution of Ax ~ 117 pc. Grid re- 
finement is triggered following three criteria, one for iV-body 
particles, the other two for the gaseous component. The par- 
ticle criterion requires that at least 24 particles reside in a 
cell before the refinement, so as to maintain a roughly con- 
stant number of particles per cell. The first criterion for gas 
is based on a ' quasi- lagrangian' approach according to which 
refining is triggered if the gas mass inside a cell exceeds a 
threshold of ~ 1.5 x 10 5 M p. The seco nd gas criterion is 
based on the Jeans length (|jeansl ll919T l : a cell at level £ 
is split if the local Jeans length does not exceed 4 times 
the cell dimension Ax^ e \ In the limit of our maximum res- 
olution, this choice allows us to properly follow the evolu- 
tion of features such as dynamically unstable gaseous clumps 
in the expanding ring. Gas dynamics is followed assuming 
a polytropic equation of state (EOS) with adiabatic index 



7 = 5/3, corresponding to a perfect mono-atomic gas. The 
interacting galaxies are embedded in a hot and thin medium 
with numerical density rib g — 10 



6 cm 3 , and temperature 



5 x 10° K. 



The RAMSES code includes physically-motivated sub- 
grid models for gas cooling, SF and supernova (SN) feed- 
ba ck. The code adopt s the cooling function described 
in ISutherland fc Dopital (Il993l ) for temperatures down to 
~ 10 4 K. We include a sub-grid gaseous EOS T = 
Tb(p/po) 70_1 to model the behavior of the inter-stellar 
medium (ISM), where T = 1000 K, 7o = 5/3 and po is th e 
SF density threshold (see lAgertz. Tevssier fc Moord [20111 ). 
SF is implemented in RAM SES following the prescriptions 
bv lRasera fc Tevssierl (|2006l ). The conversion of ga s density 
Pgas into star density p* is based on a Schmidt-law (|Schmidtl 
1959) of the form: 



dp, 
dt 



Pgas 
^dyn 



Pgas > PO, 



(1) 



where esf is the SF efficiency, td yn is the local dynamical 
time and po is a density threshold. We chose esf = 0.1 
and po corresponding to a numerical density no = 5 cm -3 . 
Such t hreshold den sity was selected to resolve the Jeans 
length i|Jeandll919T l for a temperature T ~ 10 3 K (i.e. con- 
sistent with the minimum gas temperat ure in our simula- 
tions ) with our best spatial resolution fsee lMapelli fc Maverl 
|2012| for details). Thus, the adopted threshold density is 
the highest that can be safely used with our spatial reso- 
lution and matches the one adopted in recent cosmological 
simulations th at reproduce realisti c Milky Way-like (MW- 
like) galaxies |Guedes et alj l201ll ) , and in simulations o f 
the formation of off-centred CRGs (IMapelli fc MaveijlibTl ) . 
Most of cosmological simulatio ns adopt a much low er density 
threshold (no = 0.1 cm -3 , e.g. lAgertz et al.1l201ll and refer- 
ences therein), that might unphysically enhance the SF (see 
e.g. th e discussion in lGuedes et al.ll201ll ). iGovernato et alj 
l|2010l) adopt no = 100 cm~ 3 (a factor of 20 higher than our 
threshold) in their cosmological simulations, but their spa- 
tial resolution is better than ours and their cooling function 
goes down to ~ 300 K. 

Our adopted value of the SF efficiency (esf = 0.1) is 
a factor of 2 — 10 higher t han adopted in recent simula- 
tions (e.g . lAgertz et al"l201ll). whi ch follow the observational 
hints by iKrumholz fc Tanl (20071) . On the other hand, we 
stress that (i) the esF parameter in our simulations is sim- 
ply phenomenological and does not have a direct relation 
to the true SF efficiency within molecular clouds (which we 
cannot s imulate realistically, see the discussion below), and 
that (ii) iGuedes et al.l l|201ll) show that the simulated SF is 

much more affected by po than by esf- 

Although the SF recipe used by IMapelli fc Maverl 

(|2012h is more elaborate than ours, we find similar results 
about the global SF history for qualitatively similar runs 
(see the discussion below). Star particles are then spawned 
following a stochastic approach l|Katz|[l"992l ). with a mass 
that is an integer multiple of mi mm ' = poAa; 3 . We stress 
that our best gr id resolution (117 pc) and our adopted 
cooling function (jSutherland fc Dopital 1199*31 ) are not suit- 
able to follow the formation of low-temperature (< 10 3 K) 
high-density (> 10 2 cm -3 ) gas clouds, whose typical scale 
is expected to be < 50 pc and that were studied in pre- 
vious simulations of multiphase interstellar medium (e.g. 
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Gerritsen fc Ickej 119971; IWada et al l 120021; iBottemal 120031; 
Tasker fc Brvanll2006l : IWada fc Normarj|2007l ; lAgertz et all 
20091 ). Therefore, our simulations do not model SF in cold 
clumps, which is physically the main channel for SF. Ac- 
counting for this would be computationally prohibitive to 
do for a wide set of galaxy interaction runs. SN feedback is 
implemented as thermal feedback due to Type II SN explo- 
sions characterized by a specific energy £ ~ 5 x 10 16 ergg -1 
ijAgertz et alj|201lh . A total energy amount E = esn rn* £ is 
deposited in the 27 cells around the star particle of mass m* 
during ~ 10 Myr after the star particle creation. We assume 
a SN efficiency esn = 0.1. 



Table 1. Adopted values of the model parameters for the target 
and intruder galaxy. These parameters do not change among runs. 



Parameter Target Intruder 



M h [Mq] 1.16 x 10 12 5.9 x 10 11 

r s [kpc] 20.0 12.0 

c 11.0 14.5 

M h [M Q ] 1.5 x 10 10 3.8 x 10 10 

a s [kpc] 0.7 0.6 

M d [Mq] 6.0 x 10 10 

R d [kpc] 3.5 

/las 0.075 



2.2 Galaxy models 

We model the target galaxy as a MW-like disc galaxy com- 
posed of three A^-body components and a gridded one: a 
dark matter (DM) halo, a stellar disc, a stellar bulge and a 
gaseous disc, res pectively. For the DM halo we c hoose the 
density profile bv lNavarro. Frenk fc White! (|l996l . NFW): 



1 Fraction of gas with respect to the total disc mass, 



Ph(r) = Pcrit <5 C 



(r/r s )(l + r/r 



(2) 



where p cr it is the critical density of the Univers^E r s is the 
halo scale radius and S c is defined as: 



5c 



200 c J 

~T~ ln(l + c) - c/(l + c) ' 



(3) 



where c is the NFW concentration, c = R2oo/r s , where -R200 
is the virial radius of the DM halo so that the virial mass 
of the NFW halo is M 20 o = 200p cr it (4/3) ty R% 00 . We as- 
sume an axisymmetric exponential profile to model both the 
gaseous and stellar discs: 

M d 



Pd(R, z) 



■ e - R/Rd sech 2 (z/z d ), 



(4) 



where Md is the total disc mass, and Rd and Zd = 0.1i?d 
are the radial and vertical scale length of both the gaseous 
and stellar discs, respectively. The gaseous disc has an initial 
constant temperature To = 15000 K and we assume that it 
is composed of a mixture of ionized hydrogen and helium 
with a mean molecular weight fi ~ 0.59. The density profile 
of the stellar bulge follows the spherical iHernquistl (|l990l ) 
model: 

M b 



Pb(r) = 



1 



(5) 



27ra| (r/a s )(l + r/a s ) 3 ' 

where Mb is the total mass and a s the radial scale length of 
the bulge. 

The intruder galaxy is a dwarf, spherical early-type 
galaxy composed of a NFW DM halo and a central stel- 
lar bulge devoid of any gaseous component. We chose not to 
simulate a gas-rich intruder because, in this paper, we are 
interested in studying the evolution of the gas of the target 
separately. 

The parameters for the initial conditions (ICs) of the 
two galaxies are shown in Table [T] ICs are generated 
by us ing the code described in IWidrow. Pvm fc Dubinskl 
(2008). The code generates self-consistent disc-bulge-halo 



1 p crit = 3g 2 /(87r G), where we assum e a Hubble constant Hq 
71 kms-'Mpc" 1 l lLarson et alj|201ll) . 



Af gas /(A/ gaa + M*). 



galaxy models derived from ex plicit distribution functions 
for each component (see also iKuiiken fc Dubinskil 1 19951 : 
IWidrow fc Dubinski|[2005l ). We use ~ 1.5 x 10 6 particles for 
the DM halo of the target galaxy, while the target disc and 
bulge are composed of ~ 3.7 x 10 5 and ~ 1 x 10 5 particles, re- 
spectively. The intruder DM halo is composed of ~ 7.8 x 10 
particles, while its stellar bulge has ~ 2.5 x 10 5 particles. 

The intruder-to-target mass ratio is ~ 1/2 for both 
the dark and baryonic components. Such a high value in- 
duce s a strong perturbation on the target disc (|Gerber et al.l 
1996). This choice is consistent with observations, indi- 
cating that the can didate intruder-to-target mass ra tio is 
generally > 0.2 (e.g. lAppleton fc Struck-Marcelllll996l ). We 
adopt a fraction of gas with respect to the total disc mass 
/g as = M gas /(Mg as + M^) — 0.075, consistent with obser- 
vations of MW-size d late-type galax ies in the local Universe 
jZhang et al.ll2009l . rEvoli et al.ll201lh. We do not e xplore the 
effects of varying / ga s fsee lMapelli fc M ayor 2012 for this is- 
sue). 



2.3 Interaction parameters 

In this paper, we investigate the effects of the interaction pa- 
rameters on the formation and evolution of CRGs. In par- 
ticular, we consider the following parameters (represented 
schematically in Fig. [TJ. 

(i) The initial relative velocity Vo is the relative velocity 
of the centre of mass (CM) of the intruder galaxy with re- 
spect to the CM of the target galaxy, at the initial distance 
between the two galaxies. 

(ii) The impact parameter b is the component perpendic- 
ular to Vo of the initial distance vector between the CMs of 
the two galaxies. 

(iii) The inclination angle ■& is the angle between the sym- 
metry axis of the target disc and Vo. 

We run a set of Af-body/AMR simulations exploring 
two different values of both & and Vo, and three different 
values of b. The initial parameters of our simulation suite 
are summarized in Table [2] We initially place the CM of 
the target galaxy at rest in the centre of the box; the disc 
lies in the xy plane. The CM of the intruder galaxy has an 
initial distance Do from the target CM equal to the target 
halo virial radius, Do = fljo™' = 220 kpc, and an initial 
velocity Vo; Vo lies in the yz plane. The time elapsed from 
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Table 2. List of performed simulations and of their interaction 
parameters. 



Label 


g 


6 


IVnl 
1 v U 1 




i j 


fkDcl 


[km s ] 


aUDUvzo 


0.0 


0.0 


250.0 


a0b0v65 


0.0 


0.0 


650.0 


a0b5v25 


0.0 


5.0 


250.0 


a0b5v65 


0.0 


5.0 


650.0 


a0bl0v25 


0.0 


10.0 


250.0 


a0bl0v65 


0.0 


10.0 


650.0 


a30b0v25 


30.0 


0.0 


250.0 


a30b0v65 


30.0 


0.0 


650.0 


a30b5v25 


30.0 


5.0 


250.0 


a30b5v65 


30.0 


5.0 


650.0 


a30bl0v25 


30.0 


10.0 


250.0 


a30b!0v65 


30.0 


10.0 


650.0 




Target disk 



Figure 1. Schematic representation of the ICs, showing Vo, b 
and 1?. See Section 2.3 for details. 



the beginning of the simulation to the galaxy collision is 
300 — 600 Myr, depending on Vo; we let the models evolve 
till ~ 200 Myr after the interaction time, i.e. the moment of 
closest approach. 

The set up of the encounters ranges from axisymmet- 
ric (6 = kpc, ■& = 0°) to highly off-centre and inclined 
(6 — 10 kpc, # = 30°) interactions. The smaller value for 
| Vo| = 250 kms -1 corresponds to a marginally bound in - 
teraction, consistent with observations (jMadore et al.f 2009). 
The higher value | Vol = 650 kms" 1 repres ents a very fast 
interaction within a massive galaxy groups l|Carlberg et al.l 
l200ll . lFaltenbacher fe Mathewdl2007h . 



3 RESULTS 
3.1 Morphology 

Fig. [2] shows the face-on view of the simulated target galax- 
ies at time r r i ng ~ 50 Myr after the interaction. Images show 
the mass-weighted gas density maps. Regardless of the in- 
teraction parameters, the ring-like features of the runs have 
a typical radius i? r i ng ~ 6 — 10 kpc (~ 2 — 3Rd). This value 
slightly decreases for fast encounters, but the effect is small 
because the ring expansion velocity is mainly influenced by 
the mass of the per turber, in agreement with the results of 
iGerber et all (|l996l ). 

Simulations with 6 = kpc show circular or nearly 
elliptical rings, depending on the value of the inclination 
angle. The gaseous component of the ring is quite broad 
(with a width Aiding ~ 1.5 — 2 kpc) and clumpy, due 
to the combined effects of cooling, SF and SN feedback. 
Inside the ring we observe radial gaseous spokes that are 
thicker and slightly less developed for the encounter with 
| Vo| = 250 kms - , and the presence of a gaseous and star- 
forming secondary ring that recently started propagating. 
Observations of the Cartwheel galaxy do not support our 
finding that spokes are mainly composed of gas: no HI or Ha 
emission is found in the spokes and in the inner ring (|Higdonl 
ll995l.lHigdonlll996h . However, recent observations of NGC 



922 (Pellerin et al. 



2010) show Ha emission and the presence 



of young star clusters along the incomplete ring as well as 
in the central bar and in correspondence to possible spokes. 
Also the inner ri ng of Arp 10 has Ha emission and seems to 
be a site of SF 1 Charmandaris. Appleton fc Marstonlll993l . 
iBizvaev et al l |2007| V Therefore, the amount of gas in the 
spokes and in the central regions of CRGs may depend on 
the pre-encounter distribution of gas (ISmith et al.ll2012r ). 

Off-centre collisions with 6 = 5 kpc a nd 6 = 10 kpc pro- 
duce r esults similar to those reported bv lMapelli fc Mayer] 
(|2012h . The asymmetry of the interaction induces the dis- 
placement of the nucleus in the direction of the impact. The 
ring is also asymmetric and composed of two main struc- 
tures originating from the displaced nucleus and connected 
at the end: a wide arc and a straighter arm. The wide arc 
is generally clumpy (e.g. run a30b5v25 or a30bl0v25) and 
resembles the gravitati onally unstable tidal tail produced b y 
galaxy interactions in IWetzsteinT Naab fc BurkertJ (|2007tl . 
For simulations with 6 = 10 kpc, the ring is loosely con- 
nected for | Vo| ~ 250 kms -1 , and not connected in the case 
of | Vo = 650 kms -1 . This trend is enhanced when the incli- 
nation angle $ is increased. Our results suggest an approxi- 
mate upper limit 6 < 3Rd to produce a connected ring from 
an off-centre high-speed collision. Another common feature 
of all the runs with 6 > kpc is the development of a stel- 
lar and gaseous bar-like feature in correspondence to the 
displaced nucleus. 



3.1.1 Disc warping 

Fig. [3] shows edge-on views of all the gaseous target discs at 
time Tving ~ 50 Myr after the interaction. Images show the 
effects of the interaction on the vertical structure of the tar- 
get disc. In all cases, the target is displaced from its original 
position in the direction of the intruder by < 35 kpc and 
< 15 kpc for slower and faster interactions, respectively. 
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Figure 2. Mass- weighted gas density maps in the xy plane for all the runs at time r r i ng ~ 50 Myr after the interaction. From the leftmost 
to the rightmost column: b = kpc, b = 5 kpc, b = 10 kpc. The top two rows show interactions with 1? = 0°, while the bottom two rows 
with •& = 30°. For each pair of rows, the first one represents encounters with |Vq| = 250 kms -1 , the second one with |Vq| = 650 kms -1 . 



Another common feature is the development of thin gaseous 
tails in the direction of the intruder passage, especially for 
inclined interactions with i9 = 30° . These tails can be some- 
times observe d as HI plumes departing from the target disc 
(|Higdonlll996h . 

Off-centre and low-speed interactions produce more 
disc vertical heating than faster and symmetric encounters. 
Moreover, the inclination angle induces a net torque on the 
target disc along the x axis. Therefore, after the interaction, 



the ring develops on a plane that is tilted with respect to the 
original plane of the disc. For the case with 6 = kpc and 
i? > 0°, we can estimate the warp analytically. Assuming 
that the intruder moves on a straight line in the yz plane 
with an inclination angle i? defined as above, only the x com- 
ponent of the torque N does not vanish on the entire target 
disc. It induces an angular momentum variation that can be 
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aDb0v25 al)b5v25 a0bl0v2.5 




y [kpc] V [kpcj y ppc] 

Figure 3. Mass-weighted gas density maps in the yz plane of all the runs at time T r ; ng ~ 50 Myr after the interaction. From the leftmost 
to the rightmost column: b = kpc, b = 5 kpc, b = 10 kpc. The top two rows show interactions with i? = 0°, while the bottom two rows 
with •& = 30°. For each pair of rows, the first one represents encounters with |Vrj| = 250 kms" , the second one with |Vq| = 650 kins" 1 . 



estimated as (see Appendix 1X1 for details): 



AJ(R,t) 




3tvG M int (Vr) 2 r 3 S(r) sin(2#) 
X 2(r 2 + a 2 + (Vr) 2 ) 5 /2 ' W 

where V is the module of the intruder velocity, a is the radial 
scale of the intruder (for which we assume a Plummer po- 
tential, see Appendix[X}, Mi nt is the intruder baryonic mass 



and £(r) is the exponential surface density of the target disc. 
The unperturbed component of the angular momentum can 
be estimated as: 

r-R 

3 {R)=e z 2-K drr 2 E(r)V circ (r), (7) 

where VcircM is the circular velocity of the unperturbed 
exponential disc. Then, we can parametrize the warping as: 
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,r, N /|AJ(ii,t)| 
VSwarp (it, t) = arctan 

V l J <)(,-«) 



(8) 



<y3warp is the average warping of the target disc at each ra- 
dial position, due to the gravitational torque exerted by the 
encounter. 

This simple treatment does not take into account the 
perturbation by the primary disc on the intruder motion 
and the propagation of the rin£|. Thus, it can only pre- 
dict the overall inclination of the disc, represented by the 
nearly constant asymptotic value that v? W ar P assumes for 
radii R > 7 kpc (~ 2i?d), as shown in the left-hand panel of 
Fig. Ul This plot shows the radial profile of </? W ar P deduced 
from Eq. ([8]) at two characteristic times (0 and 50 Myr after 
the interaction) for runs a30b0v25 and a30b0v65. For V we 
use the relative velocity between the CMs of the baryonic 
components of the two galaxies at the moment of the inter- 
action, as measured from simulations (V = 850 kms -1 for 
run a30b0v25 and V = 1050 kms -1 for run a30b0v65). We 
assume the radial scale length of the intruder is a — 0.75 kpc 
(see Appendix IA")l . For run a30b0v25, ip va .r P approaches 
the asymptotic value <£ W ar P — 16°, while for run a30b0v65 
if W!iI p ~ 12°. These values are in fair agreement with sim- 
ulations (as shown in the right-hand panel of Fig. [4j and 
represent the overall inclination of the target disc with re- 
spect to the xy plane. 

Off-centre interactions with ■& > 0° produce rings com- 
pletely warped in the vertical direction, especially for slow 
encounters (see Fig. [3|. In these cases, the loosely connected 
ring develops with a nearly helicoidal shape because of the 
coupled role played by the inclination angle, which induce 
the vertical torque discussed above, and by the asymmet- 
ric propagation of th e expanding ring, due to b > kpc 
l|Struck-Marcel]|ll99Gn . 



3.2 SF history 

We study the effects of the interaction parameters on the 
global SF histories of our simulations. In Fig. [5] we show 
the SFR as a function of time for all the runs with different 
inclination angle i9, while in Fig. [H] we compare the SF histo- 
ries for runs with different initial relative velocity Vo. Before 
interaction, the target has a low SFR ~ 0.1 — 1 MQyr -1 , 
consi stent with the o bservational Kennicutt- Schmidt rela- 
tion l|Kennicuttl ll998T l. For all models, the target experi- 
ences an enhancement of SF when the intruder crosses the 
target disc. The SFR rapidly increases from < 1 Mq yr -1 in 
the pre-encounter phase to ~ 10 Mq yr ~ during the star- 
burst, in agreement w i th observations jMavva et al.l 120051 . 
iRappaport et all l20ld . iFogartv et all 1201 ll ) Figs. 151 and © 
suggest that neither $ nor Vo influence significantly the SF 
history of a typical CRG. 

By contrast, the impact parameter b influences signifi- 
cantly the SF history in our simulations, as shown in Fig. 
[7] We systematically observe a double-peaked SF history 
for axisymmetric encounters. After the interaction, the SFR 
steeply increases, reaches a maximum and then decreases 



2 The impulse approximation adopted to derive equation ([Sj i s 
acceptable, as the total mass of the intruder far exceeds the fluc- 
tuations induced in the stellar density of the target. 



on a typical time scale tsf ~ 50 Myr. This is probably be- 
cause the ring travels rapidly across the disc during the first 
~ 100 Myr after the interaction; when it reaches the out- 
skirts of the disc, the lower gas density quenches the SF. 
The first SF bump is followed by a second episode of SF due 
to the propagation of the secondary ring. 

Instead, after the initial starburst, all off-centre en- 
counters show a nearly constant or slowly rising SFR that 
lasts more than ~ 100 Myr. The SF histories of runs with 
b — 5 kpc and |Vo| = 250 kms -1 represent a transition 
between truly axisymmetric and strongly off-centre interac- 
tions, as they still show a less pronounced double-peaked 

shape. 

These results are in agreement with IMapelli fc Maverl 

|2012j '). In particular, the SFR in their C and D models 
(with /gas ^ 0.09 and / gas ^ 0.05, respectively) reaches 
~ 20 Moyr -1 and ~ 8 Moyr -1 , respectively, during the 
first 50 Myr after the interaction. Our run a0bl0v65 is qual- 
itatively comparable to theirs and we find similar values of 
SFR ~ 8— 10 Mq yr -1 despite the different implementation 
of the SF in the adopted codes. 

None of the considered interaction parameters appre- 
ciably affects the peak of SFR. This feature is more sen- 
sitive to other paramet ers such as / gas , as pointed out by 
IMapelli fc Maverj (|2012T ). 



3.3 Synthetic photometry 

Photometric properties can be inferred from our simula- 
tions, by using synthetic stellar population models. We 
use mass-to-light ratio tabled ba sed on the iso c hrone s and 
synthetic ste ll ar po pulations of IMarigo et all |200Sh and 
iGirardi et ail |2010). These tables span the stellar age in- 
terval from 4 x 10 6 yr to 12.6 x 10 10 yr and give the mag- 
nitude in the UBVRIJHK bands for a stellar system with 
total mass 1 M q and a Chabrier initial mass function (IMF; 
IChabrier|[200lT) . We assume a constant sub-solar metallic- 
ity Z = 0.0019 = O.IZq, in agreement with observations of 
CRGs (|Fosburv fc Hawardedll977l ; iBransford et al.lll99Sft . 

In Fig. [5] the synthetic B — V colour maps obtained for 
all our simulations at time T r i n g — 50 Myr are shown. We 
calculated these maps after assigning a random age (between 
1 Gyr and 10 Gyr) to stellar particles that were in the ICs; in 
this process we do not take into account any pre- interaction 
age/metallicity gradient. For purpose of visualization, maps 
are smoothed with a Gaussian filter with a = 300 pc. 

All models show a blue ring-like feature with typical 



colour (B - V)r 



0.5 mag. The rings are sprinkled 



with bluer knots corresponding to regions of recent SF, with 
colour (B — l/)knot < 0.4 mag. The underlying disc is redder 
with an average colour (B — V)disc ~ 0.65 mag. The ring 
is characterized by a blue colour bump with respect to the 
disc, A(B — V) ~ 0.2 mag, which is consistent w i th B — V 
profiles derived by iRomano. Mavva fc Vorobvovl (|2008l l for 
a sample of 15 CRGs. We repeated the same analysis as- 
suming a metallicity Z = 0.019 = Zq, and we found similar 
results, even if with a redder (~ 0.15 — 0.2 mag) overall 
colour. 

Simulations with b = kpc show a gaseous secondary 



Sec http : // stev. oapd. inaf . it/cgi-bin/cmd. 
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Figure 4. Left-hand panel: radial profile of i/3 warp as calculated from Eq. (0 at t = Myr (the time of closest approach), and t = 50 Myr 
after the interaction for run a30b0v25 (solid lines) and for run a30b0v65 (dashed lines). Right-hand panel: gas density maps projected in 
the xz plane for runs a30b0v25 (top figure) and a30b0v65 (bottom figure). In both panel, the white line shows the theoretical asymptotic 
value of (/3 W arp- Each map is 40 X 20 kpc per edge. 



zr io 




zr io 



-150-100 -50 50 100 150 -150-100 -50 50 100 150 -150-100 -50 50 100 150 
t [Myr] t [Myr] t [Myr] 

Figure 5. Effect of the inclination angle 1? on the SF history. Each panel shows the SFR as a function of time, where t = Myr refers to 
the interaction time; the solid (dashed) line refers to the case with # = 0° (1? = 30°). The top row shows runs with | Vol = 650 krns" 1 , 
while the bottom row shows runs with |Vrj| = 250 kms -1 . From the leftmost to the rightmost column: runs with 6 = kpc, b = 5 kpc 
and b = 10 kpc, respectively. Error bars represent l-a uncertainty on the estimated SFR. 
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ring that just started to expand (see Fig. [2| and that is 
expected to form stars, but colour maps do not show an 
excess of blue colour in the inner regions of the ring. This 
is in partial agreement with observations: for example, the 
Cartwheel galaxy shows a r ed but gas-poor secondary ring 
l|Higdonlll995l . lHigdod[l99^ . In our models, the secondary 
ring starts forming stars ~ 50 Myr after, in correspondence 
to the second peak of the SFR at ~ 100 Myr after the in- 
teraction. 



Off-centre collisions show very blue knots mainly lo- 
cated in the wide arc and at the extremity of the bar-like 
nuclear region. This suggests that the asymmetry of the en- 
counter can funnel gas to the nucleus and favor SF near it , 
as observed, for example, in NGC 922 jPellerin et alj|20ld ). 



Our simulations likely underestimate the dumpiness of 
the blue knots in the ring, as we cannot resolve the giant 
molecular clouds and thus we cannot produce realistic SF 
agglomerates, such as young dense star clusters (see Section 
12.11 for a discussion about the limitations of the model). 



4 SUMMARY 

In this paper, we investigate the properties of CRGs as a 
function of the interaction parameters by means of AMR 
hydrodynamical simulations. We focus on three parameters: 
the impact parameter, the initial relative velocity and the 
inclination angle. 

All the simulations show a few common features, such 
as the radius of the ring (i? r i ng — 6 — 10 kpc) at time 
Tring ~ 50 Myr after the interaction, which is weakly influ- 
enced by the three interaction parameters described above. 
On the other hand, axisymmetric encounters produce CRGs 
noticeably different with respect to those formed by off- 
centre interactions. Axisymmetric interactions result in a 
circular primary ring followed by a smaller secondary ring, 
whereas CRGs born from off-centre collisions show asym- 
metric rings and a displaced nucleus. In the most asymmet- 
ric cases (6 = 10 kpc and t? = 30°) the ring is not connected. 
The maximum impact parameter to form a connected ring 
is b < 37?d, where Rd is the scale length of the target disc. 

Interactions with i? > 0° induce a torque on the target 
disc producing a vertical warp. In the special case of encoun- 
ters with 6 = kpc and ■& > 0° , the ring develops on a plane 
tilted with respect to the original target disc plane and we 
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show runs with 1? = 0° and # = 30° , respectively. Error bars represent 1-er uncertainty on the estimated SFR. 



propose a simple analytical treatment to predict the overall 
warping of the target disc after the interaction. 

After the interaction, the ring expansion enhances the 
SF in the target disc, reaching a ty pical SFR ~ 10 Mp yr -1 , 
in agreement with observations l|Higdon fc WallirJ 1 19971 . 
iFogartv et al J 12011ft . We show that the initial velocity and 
the inclination angle have weak effects on the peak of the 
starburst and on the evolution with time of the global 
SFR. Instead, the impact parameter has a crucial effect 
on SF: symmetric encounters cause a short-lived starburst 
episode, whereas off-centre interactions produce long-lasting 
(> 100 Myr) SF events. 

The expanding ring is the regio n mainly affected by SF , 
as suggested by obs ervations (e.g. iHigdon fc Wallinl 1 19971 . 
iRomano et alj [2008). This is confi rmed by the synth etic 
colour maps (based on the models bv lMarigo et al.ll2008l and 
iGirardi et ai1l2010ft that we compute from the simulations. 
We find an average B — V colour ~ 0.5 mag for the ring, with 
clumps of massive SF with B — V < 0.4 mag. The ring is 
typically ~ 0.1 — 0.2 mag bluer than the rest of the galaxy, 
in agreement with t he B — V colour profiles measured by 
IRomano et al.l |2008T l for a sample of 15 CRGs. 

A number of issues cannot be explained by our sim- 
ulations. For example, it is still unclear which mechanism 



regulates SF in the inner regions of the ring, since HI/Hq 
observations sug gest that seco n dary rings are gas-poor (e.g., 
the Cartwheel Irligdonl 1 19951 ; lHigdonlll996l ; but see also 
iBizvaev et all 120071 about Arp 10), while our simulations 
predict the opposite. This may be due to several reasons, 
such as the initial distribution of gas inside the disc, or the 
idealized sub-grid SF process that does not take into ac- 
count the chemical abundances and the role of Eb. Thus, fu- 
ture works will require more realistic initial gas distributions 
and/or more detailed recipes of SF (e . g. H2 regulated SF , 
iKrumholz. McKee fc Tumlinsonl |2009| ; iKuhlen et ail 12012ft 
to better constrain the SF in CRGs. Another interesting is- 
sue we did not explore (but plan to investigate in the future) 
is the effect of a gas-rich intr uder. In this r espect, we remind 
that the SPH simulations of IStruckl (|l997l ) suggest that the 
main differences between encounters with (disc-like) gas-rich 
and gas-poor intruders are related to the SF history, whereas 
the effects on the ring morphology are quite limited, even if 
some features (such as plumes) appear to be more developed 
in the gas-rich case. 
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Figure 8. B — V colour maps of all the simulations in the yz plane at time T r i ng ~ 50 Myr after the interaction. Maps are computed 
assuming a Chabrier IMF and a constant metallicity Z = O.IZq. From the leftmost to the rightmost column: 6 = kpc, 6 = 5 kpc and 
6 = 10 kpc. The top two rows show interactions with # = 0°, while the bottom two rows with $ = 30°. For each pair of rows, the first 
one represents encounters with |Vq| = 250 kms -1 , the second one with |Vo| = 650 kms -1 . 
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APPENDIX A: ANALYTICAL CALCULATION 
OF THE WARPING 

We model the intruder galaxy as a Plummer sphere with 
potential: 

GMint 1 



*tat(r) = 



Vi + W«) 2 ' 



(Al) 



where a is the scale radius of the potential and Mint is 
the intruder baryonic masfl We assume that the intruder 
moves on a straight line in the yz plane, according to 
R(t) = (0, — V t sin — V t cos $), where # is the inclina- 
tion angle and V is intruder velocity. For simplicity, we as- 
sume that V is constant in time. The target disc initially 
lies on the xy plane and has a surface density profile E(r) = 
n M it e~ T / Rd . The force exerted on an infinitesimal mass por- 

tion dm = T,(r)rdrd(j> located at x = (r cos <j>, r sin <f>, 0) of 
the target disc is: 

dF(r,<j>,t) = -dmV*tat(|x-R(t)|) = 

T,{r)rArA<i>GM iu t 



(r 2 + (Vt) 2 + 2V t r sin § sin + a 2 ) 3 / 2 

(A2) 



r cos <f> 
r sin (f> + V t sin i? 
Vt cos i? 



We can now derive the torque exerted on the same mass 
clement: 

dN(r,<M) = xAdF(r,(4,t) = 

GM int rVt E(r)rdrd(j) 



X 



where: 



C(r,t) 



2Vir sim? 
r 2 + a 2 + (Vt) 2 



(A4) 



Integrating Eq. (|A3|) over d(f>, we obtain the total torque 
exerted on an infinitesimally thin annulus of radius r. We 
point out that the y and z components are proportional to: 



2tt 



COS 



0. 



/, (1 + C{r,t) sin (f,) 3 / 2 
Therefore, the resulting torque per unit length is: 

2 - 



(A5) 



dN(r,t) 
dr 



GM int r 2 VtS(r) cos g 
' (r 2 + a 2 + 2(Vi) 2 ) 3 / 2 
.-in ado 



o (1 + G(r,t)sin0) 3 / 2 



(A6) 



Integrating this equation with the boundary condition 
N(7l = 0) = 0, we obtain the total torque: 



N(R,t) = 



, GM int r 2 ViS(r) cos 
dr— — ,.. ;„,„,„ x 



(r 2 + a 2 + 2(Vt) 2 ) 3 / 2 
2n sin^d^) 



lo (1 + G(r,t) sin 0) 3 / 2 ' 
The angular momentum J is related to the torque N by 



(A7) 



dJ{R,t) 
dt 



N(R,t). 



(A8) 



We can integrate in time Eq. (|A7|) to obtain the variation of 
J due to the torque exerted by the intruder, A J, with the 
condition linii-^-oo AJ — (i.e. we assume that for t — > -co 
R| = +oo and the torque is zero). We finally obtain: 



AJ (R,t) 



drN(R,t) = 

f R , GM int Vrr 2 S(r) cos^ 
dr / dr— — ^-t^tttt^— x 



o (r 2 + a 2 + (Vt) 2 ) 3 / 2 



(1 + C(r,r) sin <j>) 3 / 2 



(A9) 



This expression can be simplified by assuming sin$ <C 1. fn 
this case we have: 



AJ(R,t) 



R 

dr I dr x 
'o 

\2 _3 ■ 



3^GA/int(Vr) 2 r 3 S(r) sin(2i?) 

2(r 2 + a 2 + (Vr) 2 ) 5 /2 ' ^ iu J 

In this treatment, the z component of the angular mo- 
mentum J is not perturbed by the passage of the intruder. 
The unperturbed component of the angular momentum is 
given by: 



Jo (-R) =e z 2 



f R 2 

/ drr z E(r)V circ (r), 
Jo 



(AH) 



(r 2 + a 2 + 2(Vt) 2 ) 3 / 2 (1 + C(r, t) sin 0) 3 / 2 

(cosi? sin0 \ 
— cos $ cos <j) , (A3) 
sin ■& cos d> j 



where th e circular velocity Vcirc fo r an exponential disc is 
given by dBinnev fc Tremaind [20oi): 



-.{R) = lHy)Ko(y) - /i(y)^(y)] 1/a , (A12) 

where y — R/(2Rd) and Io,Ii,Ko,K~i are modified Bessel 
functions. The warping of the disc can be parametrized by 
the warping angle iy5 W ar P , defined as: 



V?war P (-R, t) — arctan 



|AJQR,t)j 
|Jo(A)| 



(A13) 



This angle represents the inclination of the disc with respect 
to the xy plane around the y axis at each radius R at time 
t due to the gravitational torque exerted by the encounter. 

This treatment of the torque induced by the intruder 
galaxy on the target disc is very approximated, as it does 
not take into account the gravitational effects of the target 
galaxy on the motion of the intruder and the propagation 
of the ring as a consequence of the interaction. Moreover, 
the intruder is modeled as a Plummer sphere for the sake of 
simplicity. However, we choose the radial scale a such that 
the mass enclosed inside a by the Plummer profile is the 
same as enclosed by the Hernquist profile: 



VW 2 



0.7495 a s 



1 



(A14) 



where a s is the scale radius of the Hernquist profile. Despite 
the approximations, this treatment roughly highlights the 
effects of a non-axisymmetric interaction with inclination 
angle ■& > 0° (see ggXT) . 



We neglect the role of the DM halos, for simplicity. 



© 2012 RAS, MNRAS OOO.riJfl3l 



